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ABSTRACT 

We report our search for variability in the Ha emission hne of the recurrent novae T CrB with time 
resolution 10-15 minutes. This is comparable with the time scale of the photometric flickering observed 
in this object. This is the first time that observations of the short time scale variation in emission lines 
have been made in this object. On two nights (990106 and 990107) we detected statistically significant 
variability (at 0.99 confidence level) in the line profile of Ha. This variability is confined in the central 
part of the emission line (±100 km s^^) although FWZI(Ha) is ^800 km s^^. The variability in 
the line profile is accompanied with variability of the total equivalent width, EW(Ha), ±8% for 
990106, and ±6% for 990107 (calculated from the mean value of EW). Assuming Keplerian motion, 
the variability is generated at a distance ~20-30Rq from the white dwarf, which is approximately the 
radius of the ring the stream of gas forms as it flows away from Li. On three other nights we are only 
able to put upper limits on the variability, AEW(Ha) ±2% for 980415, ±4% for 980802, and ±3% for 
980803. 

Subject headings: stars:individual: T CrB - binaries: symbiotic - binariesmovae, cataclysmic variables 



1. INTRODUCTION 

The short term photometric behaviour of symbiotic 
stars has been investigated in many systems (e.g. Do- 
brzycka et al., 1996; Sokoloski et al, 2001). However, 
there are no systematic searches for rapid spectral varia- 
tions. Until now the search for rapid spectral variability 
has been undertaken in only three symbiotics: CH Cyg, 
MWC 560, and RS Oph. All three stars are known to 
exhibit flickering and have coUimated outflows. Rapid 
spectral variability in the Balmer lines of CH Cyg on 
time scales of ~1 hour have been detected and is prob- 
ably connected with blobs ejected from a white dwarf, 
acting as a propeller (Tomov et al. 1996). For RS Oph, 
Sokoloski (2002) reported variability on time scales of 
hours in the HeIIA4686 line although no variability in 
H/3 was detected. No significant variability has been de- 
tected in MWC 560 (Tomov et al. 1995). 

In a few cataclysmic variables (CVs), which are closely 
related to symbiotic stars, rapid spectral changes are vis- 
ible in optical and UV lines as a result of a variable ac- 
cretion disk wind: BZ Cam (Ringwald & Naylor 1998), 
V592 Cas (Witherick et al. 2003), RW Sex (Prinja et al. 
2003), V 603 Aql (Prinja et al. 2000a), BZ Cam (Prinja 
et al. 2000b). 

The detection of variability on time scales of tens of 
minutes in symbiotic stars and cataclysmic variables pro- 
vides the motivation for this study. T CrB consists of a 
red giant and a hot component, most probably a white 
dwarf. T CrB can be classified as a cataclysmic variable, 
symbiotic star, and recurrent nova and accretes at a rate 
lxlO~**M0 yr-i (SelveUi et al. 1992). This places it 
among the cataclysmic variables with the highest accre- 
tion rates (see Warner 1995, p. 476) and for such high ac- 
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cretion rates an accretion disk wind could be expected. 
Here we report our search for signatures of variability 
in time-resolved optical spectroscopy in the Ha line of 
T CrB, with time resolution ^10-15 minutes. This reso- 
lution is comparable with the time scale of flickering (see 
for example Zamanov et al. 2004). 

2. OBSERVATIONS 

The observations were performed with the Coude spec- 
trograph of the 2.0 m RCC telescope at the Bulgarian Na- 
tional Astronomical Observatory "Rozhen" . The spectra 
have a dispersion of 0.2 A pixel"^. The normal CCD re- 
duction procedures (i.e. bias subtraction, flat flelding, 
wavelength calibration, etc) have been undertaken in the 
IRAF environment. A journal of observations is given in 
Table 1. The signal to noise ratio (SNR) achieved on the 
individual exposures is 35-65. 

The averaged spectra, normalised to the local contin- 
uum, are presented in Fig^ along with a comparison 
spectrum of the red giant HD135530. 

3. RESULTS 

3.1. Average, fractional and temporal variance spectra 

The prominent emission line in Fig. ^ is Ha. On a 
time scale of years the intensity of the Ha line varies from 
EW=0.5 A up to 35 A (lijima 1990, Anupama & Prabhu 
1991, Zamanov & Marti' 2000, Stanishev et al. 2004). 
During the search for rapid variability we have observed 
a range of different profiles. An atlas of the different Ha 
profiles is given in Stanishev et al. (2004). In the spectra 
obtained on 15th April 1998 the line is single peaked, 
and the central absorption dip is not visible. During 
August 1998, the profile is a little unusual, with 3 peaks 
or 2 absorption dips. During January 1999, the profile 
exhibits double peaks with a well defined central dip, 
which is partly due to the Ha absorption line of the red 
giant. 
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TABLE 1 

Ha OBSERVATIONS OF T CrB. In THE TABLE ARE GIVEN THE DATE OF OBSERVATIONS, THE BEGINNING AND THE END OF SPECTROSCOPIC 
OBSERVATIONS AND TRUNCATED JULIAN DAY. ORBITAL PHASE IS CALCULATED USING THE EPHEMERIS OF FeKEL ET AL. (2000). THE 
FOLLOWING COLUMNS GIVE WAVELENGTH COVERAGE, NUMBER OF EXPOSURES IN EACH NIGHT AND THE EXPOSURE TIME OF THE INDIVIDUAL 
EXPOSURES. For the EW(H«) we give the mean, STANDARD DEVIATION, PLUS THE SMALLEST AND THE LARGEST VALUES FOR EACH NIGHT. 



Date 


UTstart 


UTend 


TJD 


phase 


wavelength 


N X 


Exp. time 


EW(Ha) 




yymmdd 


hh:mm - 


hh:mm 


mid 




coverage 

[A] 






[A] 


[A] 


980415 


22:37 - 


01:43 


50919.512 


0.187 


6510 - 6620 


10 


X 15 min 


20.84 ± 0.31 


20.39 - 21.36 


980802 


18:50 - 


20:58 


51028.328 


0.665 


6510 - 6620 


12 


X 10 min 


13.39 ± 0.23 


12.72 - 13.81 


980803 


18:23 - 


21:03 


51029.305 


0.669 


6510 - 6620 


15 


X 10 min 


13.97 ± 0.25 


13.63 - 14.45 


990106 


02:59 - 


04:13 


51185.652 


0.356 


6500 - 6700 


7 


X 10 min 


9.77 ± 0.62 


8.90 - 10.51 


990107 


01:50 - 


02:51 


51186.598 


0.360 


6507 - 6707 


3 


X 20 min 


7.92 ± 0.45 


7.45 - 8.36 



1 1 1 1 1 1 


1 1 1 1 1 1 1 ill 
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Fig. 1. — The nightly average spectra for all five nights. A spec- 
trum of the red giant HD135530 (M2III) is shown for comparison 
at the bottom. 



Fig. 2. — Measurements of the EW{Ha). The vertical error 
bars indicate the errors in EW measurements, the horizontal the 
exposure time of the spectra. In the first three nights the variability 
is of the order of the estimated errors. On 990106 and 990106 it is 
more or about 5 times the error of the measurements. 



The total EW(Hq!) was measured on each exposure in 
our sample. The results are presented in Fig|21 In the 
first three nights the variability is less than three times 
the estimated errors. On 990106 and 990107 it is about 
5 times the error of the individual measurements. 

The measurements of the EW depend strongly on the 
continuum placement. The continuum was defined as 
a stright line or spline fitting, using parts of the spec- 
trum located on both sides of Ha. We used different 
wavelength intervals to produce different normalisations. 
This leads to changes of the absolute value of EW, but 
keeps the relative distributions of the points the same. 
The errors of the EW have been calculated from four 
measurements, with different normalizations. 

In order to search for the variability in the line profile 
of Ha line during an individual night, we calculate the 
fractional variance spectrum based on all the data in each 



night: 



N 



(1) 



where N is the number of spectra, fi\ is the flux at wave- 
length A of the i*^ spectrum and f\ is the flux at that 
wavelength averaged over all spectra. The a in the con- 
tinuum should correspond to (SNR)^^, and hence a peak 
in a should determine the significance of a variation. 

One other statistical technique which is particularly 
useful to accurately determine the spectral location of 
variability is the temporal variance spectrum (Prinja et 
al. 2003; FuUerton, Gies & Bolton 1996): 



{TVS)x = 



N 



{N 



{fix — fwxY 



(2) 
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where f^x is the weighted mean of the normahsed inten- 
sity, Si is the inverse of the SNR of spectrum i measured 
in the continuum, and Sq — 

The resuhs of calculations are presented, together with 
all spectra, in Fig|2| 

As can be seen, the fractional variance spectra do 
not indicate variability in the first three nights (980415, 
980802, and 980803), however the TVS indicates some 
variability in 980415. A clear peak in both a and the 
TVS is visible on 990106 and 990107 at about velocity 
«0 km s-^ 

Following Fullerton, Gies & Bolton (1996), if the TVS 
for a spectral feature is above a specified level of signif- 
icance then the null hypothesis of "no variability" can 
be rejected. The statistical distribution of the TVS 
is governed by the the reduced chi-squared distribu- 
tion with N-1 degrees of freedom and scaled with SNR 
(TVS^ SqXn-i)- For the nights considered here, the 
statistical significance of 1% corresponds to TVS^/^=3.6- 
5.0%, depending on the number of spectra and SNR. In 
the central part of the line the TVS is well above this 
level for 990106 and 990107 (FigEJ). The peaks in a 
correspond practically to the same confidence level of 3 
standard deviations from the mean. This clearly shows 
that in these two nights there is statistically significant 
variability in the Ha line profile, at a confidence level of 
0.99, at a time resolution of ^10 min. 

3.2. Variability in January 1999 

For the two nights with detected variability, we sub- 
tract from the individual spectra the average spectrum, 
the minimum spectrum and the maximum spectrum. 
The average spectrum is defined as the average of the 
normalised flux for every pixel. The minimum (max- 
imum) spectrum is defined as the minimal (maximal) 
value for each pixel. We do not know a priori whether 
we have variability due to a pure emission, pure absorp- 
tion, or an admixture of both. The results for 990106 
are presented in Fig. 0] As can be seen, variability ex- 
ists in all three panels. This indicates additional absorp- 
tion/emission with equivalent width of approximately 1 
A, located in the central part of the line. 

To exclude the possibility that the observed variability 
is an artifact of data reduction we performed tests with 
data processing and normalisation, which showed that 
the detected variability on 990106 and 990107 is not an 
artificial result. Figl^b shows the results from two differ- 
ent extractions and normalisations which have identical 
a peaks. 

Although the sky conditions looked like photometric 
during the time of the observations we can not exclude 
the possibility that the variable atmospheric lines and 
small wavelength shifts from spectrum to spectrum can 
result in spurious variability. To test the variability these 
effects could produce, we artificially introduced random 
shifts in wavelength up to 3 km s^^[the biggest shifts we 
detected were <0.25 px (2.5 km s^^)]. We also isolated 
athmosperic absorption lines from a high SNR spectrum 
of Spica {a Vir) and introduced variable atmospheric ab- 
sorption in the spectra. We find that both efects can 
produce variability, but this variability is lower than the 
noise in TVS and a and can not account for the observed 
peaks. 



It is known (i.e. Zamanov & Bruch 1998) that the hot 
continuum, which fiickers on time scales of minutes with 
amplitude 0.2-0.3 mag in U, contributes about 15% of the 
average flux in the V band. To check whether the vari- 
ability detected in Ha is due to the changing continuum 
level caused by the flickering of the hot component, we 
simulate a in such a case by supposing that the hot con- 
tinuum flickers with amplitude ±25-30% and contributes 
15% of the flux in the continuum around Ha. The am- 
plitude is chosen to be slightly higher than the highest 
one observed in the U band after the subtraction of the 
red giant contribution (see Zamanovfe Bruch 1998). To 
this variable hot continuum we add constant red giant 
and Ha emission. The simulated a and TVS are shown 
in FigO The simulated quantities are different from the 
"observed" ones as (i) they have smaller amplitude and 
(ii) they are wider and spread over the whole FWZI of 
the Ha emission, while the observed ones are located in 
the central part of the line only. This indicates that de- 
tected variability is not due to the continuum variability 
but is intrinsic to the line. 

4. DISCUSSION 

The FWZI of the Ha emission line is 800 km s^^, with 
wings of the line formed in the immediate vicinity of the 
white dwarf (see also Stanishev et al. 2004). However, 
the detected variability appears in the central part of 
the line and is confined to /SV « ±100-120 km s~^ from 
the systemic velocity. This indicates that this variability 
originates in the outer parts of the Ha emitting region. 
Assuming Keplerian motion, a 1.4 Mq white dwarf, and 
inclination 67" (i.e. Stanishev et al. 2004 and references 
therein) this corresponds to a distance >2O-25R0. 

Mass accretion in T CrB is occurring via Roche lobe 
overfiow. As the stream of gas fiows away from Li it 
forms a ring (see Verbunt 1982). The position of the 
ring approximately defines the outer edge of the accre- 
tion disk. The radius of the ring can be estimated from 
(Hessman & Hopp 1990): rja = 0.0859(g)-°-426, where 
a is the semi-major axis and q is the mass ratio of giant 
to white dwarf. For T CrB q=0.82±0.10 and a « 210 
Rq (Stanishev et al. 2004), which places the ring at 20 
R0. At this location, the Keplerian velocity is about 115 
km s^^, which is similar to ISV of the detected variabil- 
ity. However, the Keplerian time scale is about 18 days, 
which is not comparable with the timing of our observa- 
tions. This suggests that the variability is not produced 
by an inhomogeneous disc structure. 

The observed behaviour of Ha emission from T CrB 
with time resolution 10-15 minutes is different from that 
of CH Cyg and of the variable winds observed in some 
CVs. We do not see ejected blobs, like those observed 
in CH Cyg (by Tomov et al. 1996), nor high velocity 
components observed in the variable disk winds of some 
CVs (Prinja et al. 2003). It might be that they are 
not detectable in our Ha profiles, even though they may 
exist and could be detected in the UV, where the hot 
component is the dominant source of radiation. 

One possible explanation for the detected variability 
could be variable absorption from the wind of the giant 
or from the accretion stream. It is worth noting, that no 
variability in the line is observed at orbital phases where 
the line of sight passes through the expected position of 
maximum absorption from the red giant wind. Another 
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Fig. 3. — The Ha profiles of T CrB for 5 nights of observations, normalised to the local continuum. The X-axis is in heliocentric velocity. 

a and the TVS for all nights are plotted. The dotted horizontal line indicates a 1% threshold for the detection of variability. For first 3 
nights a docs not indicate any variability, while the TVS is marginally above the 1% significance level on the night of 980415. On 990106 
and 990107 well defined peaks, both in a and in the TVS, are clearly visible in the central part of the line, corresponding to velocity -100 
to +100 km/s. This demonstrates that there is a statistically significant variability in the profile of Ha line with time resolution ~10 min. 
The probability that this is a real variability is >0.99. 



cause could be additional emission from the area where 
the accretion stream hits the disk (the hot spot). It is 
possible that a variable component with EW<0.5A al- 
ways exists, but when the lino is strong it is lost in the 
line (e.g. 980415). 

An additional plausible explanation is that the variable 
continuum arising from the central part of the accretion 
disk (where the flickering in U arises) changes the ionisa- 
tion state of the outer regions of the (flared) disk, from 
which the central part of the emission lines arise. This 
possibility deserves further exploration. 



5. CONCLUSIONS 

We searched for rapid Ha variability in the spectra of 
the recurrent nova T CrB. On 2 nights out of 5, we detect 
statistically significant variability at time resolution of 
10-20 minutes. The detected variability is confined to 
the central part of the line profile at ±100 km s""'^, which 
suggests that variations are produced in the outer parts 
of the accretion disk (>20 R©). On these two nights 
there is evidence for variability of the total EW(Ha) (in 
all cases it is calculated relative to the average value, see 
Table 1) with amplitude ±8% for 990106 and ±6% for 
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Fig. 4. — Night 990106: from tlie individual spectra have been 
subtracted tlie average spectrum (left panel), the minimum (mid- 
dle panel), and the maximum spectrum (right panel). Zero of 
the velocity axis corresponds to the systemic velocity 7=-27.79 
km s-l(Fekel et al 2000). 
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Fig. 5. — (a) Ha profiles observed on 990106 (the same as 
in Fig|3J. (b) cr for 990106, solid and dashed lines are for two 
deliberately different extraction and normalisation procedures, (c) 
simulated cr as a result of hot continuum variability, (d) simulated 
TVS as a result of hot continuum variability. During the simulation 
we assumed that the hot continuum contributes 15% of the flux on 
average and it varies with ±30% (dashed line) and ±25% (solid 
line). Zero of the velocity axis corresponds to the systemic velocity 
7- 
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990107, which is approximately five times the error of 
the individual measurements. 

For the other 3 nights of observations we do not detect 
line profile changes. We place upper limits on the vari- 
ability of the total EW(Ha) AEW(Hq) ±2% for 980415, 
±4% for 980802 and ±3% for 980803. 

More extensive spectral observations with better time 
resolution, combined with photometry are necessary to 
define what influences the Ha variability and whether 
this variability is directly connected with the flickering. 
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